Bose-Einstein condensation in an optical lattice 
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In this paper we develop an analytic expression for the critical temperature for a gas of ideal bosons in a 
combined harmonic lattice potential, relevant to current experiments using optical lattices. We give corrections 
to the critical temperature arising from effective mass modifications of the low energy spectrum, finite size 
effects and excited band states. We compute the critical temperature using numerical methods and compare to 
our analytic result. We study condensation in an optical lattice over a wide parameter regime and demonstrate 
that the critical temperature can be increased or reduced relative to the purely harmonic case by adjusting the 
harmonic trap frequency. We show that a simple numerical procedure based on a piecewise analytic density of 
states provides an accurate prediction for the critical temperature. 
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I. INTRODUCTION 

Bosonic atoms confined in optical lattices have proven 
to be a versatile system for exploring a range of physics 
Q H S BH II 0] , exemplified by the superfluid to Mott- 
insulator transition S(J. In the superfluid limit a condensate 
exists in the system and experiments have explored its proper- 
ties, such coherence OaSllfllllll, collective modes 1 12], and 
transport |0, [H, [3l • While several experiments have con- 
sidered the interplay of the condensate and thermal cloud at 
finite temperatures ifTol [T2I1 . the nature of the condensation 
transition itself remains to be examined. 

For the 3D Bose gas significant theoretical attention has 
been given to the condensation transition. The ideal uni- 
form gas has the well-known critical temperature Tco ^ 
(N/Vy^^, although it is only recently that consensus has 
been reached on how s-wave interactions shift this result 
El [M [H [H El El- In the experimentally relevant har- 
monically trapped case the ideal transition temperature scales 
as Tco ^ ujN^/^ in the thermodynamic limit. Finite size 1I21I1 
and interaction effects (at the meanfield level) ||22|] give im- 
portant corrections, and including them is necessary to obtain 
good agreement with experiment ll23ll (also see ll24ll25ll ). 

While the occuiTence of condensation in a lattice is hardly 
surprising (when interactions are small), there are few the- 
oretical predictions for the condensation temperature or be- 
haviour. For the idealized case of a (uniform) translationally 
invariant lattice Kleinert et al. [26] have made predictions that 
a re-entrant phase transition will be observed with varying in- 
teraction strength. Ramakumar ef aZ. |27] have also examined 
interaction effects in the translationally invariant lattice, and 
have explored the critical temperature dependence on lattice 
geometry. 

In experimentally produced optical lattices the periodic po- 
tential is always accompanied by a harmonic potential, pro- 
duced by the focused light fields used to make the lattice and 
sometimes enhanced by magnetic trapping (e.g. see ll28ll ). We 
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Figure 1 : Schematic diagram showing the combined harmonic lattice 
potential considered in this paper. 



refer to this experimentally realistic potential as the combined 
harmonic lattice potential (see Fig. [Hi. We are aware of two 
numerical studies that have considered finite temperature con- 
densation in the combined potential ll29l[30ll . Wild et al. ll29ll 
have considered a quasi- ID system and examined the effect 
of interactions on the transition temperature using a meanfield 
approach. Ramakumar et al. fid\ used numerical studies to 
examine condensation and thermal properties for the ideal gas 
limit. All of these studies ll26l l27l [29ll30tl have used a tight- 
binding description (or Bose-Hubbard model) that neglects 
the role of higher vibrational bands, and can only be applied 
when the lattice is sufficiently deep and the atoms are suffi- 
ciently cold. Going beyond the tight-binding approximation 
Zobay et al. Isill have used meanfield and renormalization 
treatments to consider the effects of interactions in a uniform 
system with a weak one-dimensional translationally invariant 
lattice (depth less than a recoil energy). 

We also note several studies showing how adiabatic varia- 
tions of the lattice depth might be used to prepare a condensate 
or reversibly condense the system ll32l O and recent papers 
debating the use of interference peaks as a signature of con- 
densation [28, 34, 35J. 

The central difficulty in calculating the properties of a quan- 
tum gas in the combined potential is that the spectrum has a 
rich and complex structure. Several articles have considered 
aspects of this system ll36l IstI [38l I39I1 for the case of one or 
two spatial dimensions. The first study we are aware of is a 
tight-binding description of ultra-cold bosons by Polkovnikov 
et al. [40]. Refs. jM Q have made detailed studies of 
the combined potential spectrum (also within a tight-binding 
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description), and more recently closed-form solutions were 
given by Rey et al. Il39i1 . In Refs. iIstI |4i1i an ideal gas of 
fermions in a ID combined potential was examined without 
making the tight-binding approximation. All of these stud- 
ies have confirmed that, for appropriate parameter regimes, 
parts of the single particle spectrum will contain localized 
states. This is in contrast to the translationally invariant sys- 
tem, where inter-atomic interactions or disorder are needed 
for localization to occur (e.g. see ll9l l42ll ). Experiments with 
ultra-cold (though non-condensed) bosons ll43ll have provided 
evidence for these localized states. 

In this paper we present a theory for an ideal Bose gas in a 
three dimensional combined harmonic lattice potential. Our 
treatment includes excited bands, and is thus valid at high 
temperatures where tight-binding descriptions fail. Central to 
our approach is the division of the spectrum into two regions: 
(1) A low energy region consisting of extended oscillator- like 
states, modified from those of the harmonic potential by the 
low energy effective mass. (2) A high energy region contain- 
ing localized states that have an energy related to the local po- 
tential energy, and also includes modes in the first vibrational 
excited bands. The division between these regions is made by 
the use of an effective Debye energy for the system. 

In current experiments with optical lattices inter-particle in- 
teractions are typically important, at least in determining the 
near zero temperature manybody ground state. Interaction ef- 
fects near the critical region have yet to be examined for the 
combined potential (although the quasi- ID case is examined 
in 1 29]), and our ideal gas results will be a useful basis for 
comparison with future studies. 

We begin in Sectionllllbv introducing relevant energy scales 
and describe an approximate analytic spectrum and density of 
states in the combined potential. In Sec. |III]we compare our 
analytic density of states to the results of full numerical calcu- 
lations to justify the validity regime of our analytic approach. 
We then derive an analytic approximation to the critical tem- 
perature in the combined potential and calculate corrections 
resulting from the low energy spectrum and the effects of ex- 
cited bands. Those results are compared to full numerical 
calculations to assess their accuracy and validity. Finally in 
Sec. HV] we present some general numerical results for the 
condensation phase diagram in the combined potential. We 
show how varying the harmonic confinement can be used to 
raise or lower the transition temperature relative to the pure 
harmonically trapped case. 
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Figure 2: Schematic diagram of important energy scales, (a) Band 
structure of a translationally invariant ID lattice of depth AEr. En- 
ergy scales eo (black horizontal line) , wo (grey dashed line) and ei 
(dash-dot line) identified (see text), (b) Correspondence of these en- 
ergy scales to the combined potential. Small thick horizontal lines 
indicate the energies of ground band localized states. 



where h is the reciprocal lattice vector, and {Vi, V2, T^} are 
the lattice depths in each direction. It is conventional to de- 
fine the recoil energy Er — hiOR = /8ma^ as an energy 
scale for specifying the lattice depth, where a — 27r/6 is the 
direct lattice vector. Properties of the single particle spec- 
trum have been discussed by several authors, e.g. see Refs. 
iS, m El m B. Here we use the results of these stud- 
ies to suggest an approximate (piecewise) analytic density of 
states for the combined lattice appropriate for determining the 
critical temperature. 

We begin in the next subsection by defining a set of useful 
quantities that will be crucial for developing approximations 
to the spectrum in different regimes. These quantities can be 
determined from solutions of the much simpler translationally 
invariant lattice (i.e. Eq. ([T) with all ujj — 0) or from analytic 
approximations valid in the tight-binding regime. 



B. Energy scales from the translationally invariant lattice 

A one-dimensional depiction of the important energy scales 
is given in Fig. |2] There we show the one dimensional band 
structure [Fig. |2ja)], and indicate several energies that we dis- 
cuss further below. 



II. FORMALISM 



1. Block state parameters 



A. Single particle Hamiltonian 

We consider the case of a single particle Hamiltonian of the 
form 



H = 



3 



2m ^ 



sin' ( ^) + 



(1) 



The quantity e„ refers to the minimum (Bloch state) en- 
ergy of the band with n- vibrational quanta, and in the full 3- 
dimensional case we will use the notation e„j to denote the 
particular n-th excited band by specifying additional quan- 
tum number(s) j. Here we will only refer to a few of these 
band minimum energies: eo is the ground state energy of the 
translationally invariant lattice and gives a lower bound for the 
ground state energy when the harmonic trap is added; eij is 
the lowest energy of the first vibrational excited state, with the 
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quantum number j = 1,2,3 used indicate that the vibrational 
excitation is directed along the -direction; We use 62 to in- 
dicate the energy above which higher excited bands become 
accessible IWal. 



C. Spectrum and density of states in the combined harmonic 
lattice potential 

1. Low energy spectrum (e < -Ble) 



2. Wannier state parameters 

We define wq as the energy of a localized Wannier state in 
the ground band. Wannier states are defined as a Fourier trans- 
form of the ground band Bloch states (e.g. see '45']), and 
as such its energy is the mean energy of all the ground band 
Bloch states (see Fig. HJa)). The tunneling between neigh- 
bouring Wannier states in the -direction is characterized by 
the tunneling matrix element Jj . This is given by the Fourier 
transform of the ground band Bloch dispersion relation along 
direction Xo. 



3. Effective mass 

Intermediate between the extended Bloch states and local- 
ized Wannier states we will need to describe finite extent 
wavepackets in the ground band. A convenient quantity for 
doing this is the effective mass at zero quasimomentum, m*, 
defined as 



1 



(2) 



q=0 



where eo(q) is the dispersion relation of the ground band and 
q is the quasimomentum. We note that the effective mass may 
be different along each direction. 



4. Tight-binding expressions 

All of the above quantities are easily obtained from calcu- 
lations of the translationally invariant lattice or equivalently 
from the well-known properties of the Matthieu functions. 
However the tight-binding limit, which should be applica- 
ble when Vj > hEji, yields several simple analytic expres 



sions for these quantities. In the appendix of Ref. [^^ 
an approximation for the energies of the excited bands is 
developed using a harmonic oscillator approximation. Us- 
ing those results we obtain eo ~ ^/^^j^R^ ^ jEji, 

eij « {2y^VjER - Er) + eo, and 63 w 4:y/mm{V.j}ER - 
+ eo. The tunneling matrix element can also be cal- 
culated using the harmonic oscillator approximation, giv- 
ing J^- « ii/^){V,/ER)^/^eM~2y/Vi/E]^,)ER. In the 
tight-binding limit the ground band dispersion relation is ap- 
proximately given by eo(q) ~ 4Jj- sin^((7ja/2) + eo, 
where q is the quasimomentum. From this we obtain expres- 
sions for the Wannier energy wq ~ cq + J^j ^•/j, and the 
effective mass m* — h^/2Jja^. 



The low energy states in the lattice are extended wavepack- 
ets, with a harmonic oscillator envelope. Indeed, the spectrum 
is that of a harmonic oscillator but with the frequency modi- 
fied by the effective mass, ^^j^ yjm/m*ujj ll39ll 



I.e. 



eLE(n) = eo 



(3) 



where the {rij} are non-negative integers. This low energy 
description is valid for quantum numbers in the range < 



Uj < Afj where 



(4) 



(see II39II ) as for values of rij greater than Afj the states become 
localized (see below). 

The density of states for these modes is given by 



ffLE(e) = 



2^3^^ 



eo < e< E 



LE, 



(5) 



where uj* ~ ^/u 



is the geometric mean of the effec- 



tive trap frequencies. The boundary of the rectangular region 
of {nj}-space, where the low energy description is valid, does 
not correspond to a well-defined energy cutoff. We introduce 
an effective Debye energy, £^le, such that a total of JV1M2M3 
low energy states would he below this energy. A simple cal- 
culation yields 



E, 



LE 



4v^ft 



J 



Co, 



(6) 



where J and m* are the geometric means of the tunneling 
matrix elements (Jj) and effective masses (iri*) respectively. 
Since E'le depends on J it is exponentially suppressed to- 
wards eo as the lattice depth increases. 

The ground state energy of the combined potential, corre- 
sponding to the state in which the condensate forms, is given 
by Eq. ((Sji with ?ii = n-z = 713 = 0, i.e. 



eo 



(7) 



We see that the effect of the harmonic confinement is to shift 
the ground state energy upward from that of the translationally 
invariant lattice i.e, eg. However, eg will still be less than wq 
if the harmonic potential is less confining than a single lattice 
site. 



2. Localized spectrum (e > EhE) 

The next part of the spectrum consists of localized states, 
arising because the offset in potential energy between lattice 
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sites near the classical turning point exceeds the respective 
tunneling matrix element. The nature of these states and the 
derivation of their respective density of states is treated fully 
in Ref. [44], but we briefly summarize those results here. 

The energies of the locahzed states are given by the local 
potential energy 

eLo(n) = ]^'ma'^{ujln\ + uolnl + i^lnl) + Wq, (8) 

where {rij} are (positive and negative) integers that specify 
the site where the state is localized. As these states localize 
to approximately a single lattice site, their energy offset from 
the lattice site minimum (i.e. ^ma? J^j ^j'^^f) is given by the 
Wannier energy wq. Schematically these states are indicated 
in Fig. |2lb) as horizontal rungs in each lattice site (recalling 
that for eLo(n) < E'le tunneling delocalizes these states). 

This description is valid for all energies above i?LE, how- 
ever for sufficiently high energy scales additional vibrational 
states become available. Here we will also approximate these 
excited band states using a localized description, i.e. 

(9) 

where we have approximated the zero point energy of these 
states as eij. Note that because the vibrational excitation may 
be directed along any coordinate direction we have three first 
excited bands to include. 

The density of states for the spectra given in Eqs. (O and 
©is 

3 

5'Loc(e) = go(e - wq) + ^go(e - ey), ^^le < e < 62, 

(10) 

where 

^/o(e) = i|(g)''V^(.), (11) 

u) — -^ujiUJ2i^3, (e.g. see ll44l l46ll ). and 6{e) is the unit 
step function. We also note that the case of a general (non- 
separable lattice) has the same density of states if we instead 

identify En ~ /8mVc^^ , where Vc = \sli- (a2 x aa)! is the 
unit cell volume and {ai, a2, as} are the direct lattice vectors. 

The localized states description of the first excited band is 
the most severe approximation we make for the combined po- 
tential spectrum, particularly because the lowest energy states 
of the excited bands will also be harmonic oscillator-like. For 
deep lattices the tunneling rates for the ground and excited 
bands are small and the localized description improves. For 
the theory we develop here, the first excited bands are as- 
sumed to be a rather large energy scale compared to the criti- 
cal temperature and this approximation should be adequate. 

3. Bare oscillator states (e > -Ebarc) 

At sufficiently high energy scales the lattice has only 
a small effect on the energy eigenstates and the spectrum 



crosses over to bare oscillator states. This cross-over oc- 
curs when the single particle energies exceed the lattice depth 
which in 3D we can take as the sum of the lattice coefficients 
^'baio = J2j ^j- The bare oscillator spectrum is of the form 
given in Eq. (|3]i but with the bare trap frequencies, i.e. 

ebaro(n) = ev + (j^J + ^ ' ^^2) 

where nj are non-negative integers. The constant ey = 
i J2^=i ^ is spatial average of the lattice potential and 
gives the shift of the high energy spectrum. The density of 
states is given by 

Share (e) = ^.nL , £ > £^bare. (13) 

For the parameter regimes of interest (lattices with depths 
greater than a few recoils) i?barc is sufficiently large that the 
bare oscillator states do not play an important role in deter- 
mining the condensation properties for the system. 



4. Intermediate energy region 

In sufficiently deep lattices many excited bands may be 
bound by the lattice, and will contribute to the density of 
states. In this case for energies greater than 62 and less than 
Ehaxc, the various density of states we have already outlined 
above will be inadequate. It is difficult to provide a reliable 
analytic description of these excited band contributions for 
several reasons: (1) Anharmonic effects of the lattice make 
predicting the locations (i.e. e„j) of these bands difficult. (2) 
The tunneling between sites in excited bands is much larger 
and worsens the localized state approximation. This necessi- 
tates an effective mass modified harmonic oscillator treatment 
(c.f. Eq. ^) that crosses over to localized states at higher en- 
ergies. Furthermore, large asymmetry between directions can 
occur depending on the orientation of the vibrational excita- 
tions of each band, making any form of Debye approximation 
of limited use. 

Here we do not treat these higher bands analytically. For 
typical experimental parameters the energy scale of these 
modes (i.e. 62) is well above kTc, and a complete descrip- 
tion is not required 

D. Full numerical solution 

To test the predictions of this paper we have made a full 
numerical solution for the single particle eigenstates of Eq. 
([TJ. To do this we use the separability of the Hamiltonian 
to convert this eigenvalue problem to a set of three ID prob- 
lems. Because the harmonic potential is quite weak (typically 
LOj ^ 0.01 — O.OStJfl' in experiments), a large number of lattice 
sites need to be represented to find eigenstates up to a conve- 
nient maximum energy (usually E'max ~ Eg + 25i?fl), chosen 
so that the density of states we construct will be useful for 
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Figure 3: Comparison of numerical smoothed density of states g{e) 
(dots) to analytic density of states g(e) (black lines) for a 3D com- 
bined harmonic lattice potential. The lattice depth parameters are the 
same for each direction, i.e. Vj = Vq for j = 1, 2, 3. For reference, 
the characteristic energy scales _Ele , eij, 62 and -Ebarc are shown 
as gray curves. Isotropic harmonic trap taken with u = O.OltJfl. All 
energies are measured relative to tg. 



ffLE(e) 
5Loc(e) 



eo < e < Sle, 
Ei^E < e. 



(15) 



Of course this result can only be expected to furnish a good 
description for e < 62- In the context of current experiments 
this energy range should be sufficiently large that this piece- 
wise density of states will be useful over a broad parameter 
regime. E.g. for *^Rb in a lOEj^ deep lattice (a — 425nm), 
we have that (e2 — eo) ~ lOEn ^ k x 1.5^K. 

We make a few observations regarding the results in Fig. |3] 

1 . For shallow lattices i?barc may be sufficiently small that 
the transition to bare oscillator states occurs before £'le 
is reached. Indeed, the density of states is mostly har- 
monic oscillator like (i.e. g ~ e^) for lattice depths less 
than AEji, and for this reason the analytic result is only 
shown for depths greater than this. 

2. For lattice depths less than lOEji the onset of the lo- 
calized excited band states in g(e) for e ~ eij is too 
rapid compared to the numerical results, arising because 
the lowest energy states in the excited band are har- 
monic oscillator like. However, agreement is observed 
to improve with increasing lattice depth, such that for 
Vo ^ 20-Eij the numerical and analytic results are al- 
most indistinguishable. 

3. The energy scale 62 increases quite rapidly with lattice 
depth, justifying our neglect of additional excited bands 
in the analytic density of states. 



temperatures up to about T ^ bEn/k. We use a planewave 
decomposition to represent the eigenstates of the combined 
potential, chosen because it provides an efficient representa- 
tion of the rapidly varying lattice potential. Typically of order 
10'' — 10^ planewave modes are used to represent the several 
thousand eigenstates in the energy range of interest. 

For the purposes of comparison to our analytic results, it is 
useful to construct a smoothed density of states, defined as 



dt 



ijk 



die 



.(1) 



,(2) , J3)- 



), (14) 



that gives an average number of eigenstates per unit energy 

with energies lying within Ae of e, where {e'"'"'} are the (ID) 
single particle energies (i = 0, 1, . . .) in the Xj- direction ob- 
tained form the numerical diagonalization. 



III. RESULTS 
A. Density of states 

Here we investigate the accuracy and applicability of our 
combined density of states Q and ( fTOl l by comparison with 
the smoothed density of states obtained from the full numeri- 
cal solution (see Fig. O. For definiteness, the analytic density 
of states we use is constructed piecewise from results (|5]l and 



B. Analytic prediction for the critical temperature 

As is apparent from Fig. [3] for lattices with Vj > AEn, the 
majority of the low energy spectrum is well described by the 
first term of the localized density of states ( fTOl i, and thus we 
use this term to estimate the critical temperature. 

The total number of particles in the ground band localized 
states, as a function of inverse temperature {(3=1/ kT ) and 
chemical potential (/i), is given by 

50 (e - Wo) 



de 



16 



3 

^2 



1' 



(16) 



I(J7) 



where gs{z) = /k^ is the polylogarithm function. 

Following the usual procedure lMTII we identify the critical 
temperature for condensation by taking the gas to be saturated 
ijjL wq) and setting A^loc(/3cOi wq) = N (the total number 
of atoms), giving 

0.4141 /hw^\ ./o 

Tco - (18) 

k \ujR J 

where /3oc — l/fcTco and we have used that 53/2(1) = 
C(3/2) « 2.612, with ({s) = J2t=i 1/^' *e Reimann zeta 
function. This expression has the same N^^^ dependence as 
the critical temperature for the uniform Bose gas. 
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C. Corrections to analytic critical temperature and excited band states (fTOl i are in general significant. We 

now consider the effect of these on Tco under the assumptions 
Expression (HI for Tco is based solely on the localized that {Ele - eo)/kT^o < 1 and (ea - eo)/kTco > 1. 
ground band states. The effect of the low energy states Q 

I 



1. Low energy correction 



The first correction we consider is to account for the low energy spectrum, described in Sec lIIC l1 To do this we replace 
go{e — wq) in ( fT6l ) for e < Ei,e by the low energy density of states (|5]). This changes Ni^ocif^co, l^-) by an amount 



AiVLE = 



de 



gLE(e) 



cO 



g/3co(e-eo) _ I 
2 



goje-wo) 

g/9co(«-tuo) _ 1 ' 



Wo 



■ ,.—.3 Z2 V^l 



where we have assumed that (i?LE ^ gq) <C fcTco to arrive at the last Une |l5 



(19) 
(20) 



2. Chemical potential correction 

Associated with the change in the low energy density of states is the change in ground state energy from (for the localized 
spectrum) to €g (for the low energy spectrum Q)- Replacing the saturated chemical potential by the ground state energy, i.e. 
setting — > Cg in ( [TtI i we obtain 



-32/c.«x3/2 



Wo 



irkT, 



cO 



kT, 



cO- 



(21) 



In deriving this result we have assumed that (wo — eg)/kTco <C 1, so that we can approximate the argument of the polylogarithm 
as 1 — {wo — eg)/kToc, and use the expansion 53/2(1 — x) ~ C(3/2) — l^pKx — C(l/2)a;. We note that Cl^) ~ —1-460 and the 
square root term accounts for the infinite slope of (73/2(2^) at z = 1. 

I 



3. Excited hand correction 



4. Corrected critical temperature 



As discussed in the derivation of Eq. ( fTOl i. at an energy 
scale of eij excited band states become accessible to the sys- 
tem, and contribute additional states described by the density 
of states 5o(e — eij) . The additional atoms accommodated in 
these states at Tco is given by 



Combining all the above results we arrive at a new estimate 
for the transition temperature. To do this we set 

N = Ni^ocificuwo) + AiVLE + AN^ + ANeb, (23) 

where j3ci = l/fcTd is the corrected transition temperature. 
Assuming that \Tci — Tco\ ^ Tco, we obtain 



Tci ~ Tco 



1 - -(AA^LE + AiV^ + ANeb)/N 



(24) 



AN, 



EB 



E 

7=1 ^ 



3/2 



(22) to first order in the AA^-corrections. The validity conditions 
are, as stated above, that (i?LE — eo)/kTco <C 1 and (e2 — 
eo)/fcTco ^ 1. This will ensure that all the changes (AA^) 
- (eij -Wo) /kTco are small compared to A^, however we caution that sometimes 

' due to cancellation a particular AA^ can be small even when 

the validity condition is not satisfied. 

We make the following observations on these corrections: 



AA^le: The low energy density of states tends to increase much 
where we have taken {eij — w^) ^ kT^o- Note that in cal- more slowing from its zero point than the localized den- 

culating this term we have summed over all contributing first sity of states does. Thus in replacing go{e) by .gLE(e) in 

excited bands. ( fTSI l, the number of states at low energy and hence the 
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number of atoms in the saturated thermal cloud both 
decrease. This leads to an increase in the critical tem- 
perature. 

ANfj,: The downward shift of the chemical potential when we 
change the saturated chemical potential from wq to Eg 
leads to a decrease in the number of atoms in the satu- 
rated thermal cloud, and hence an increase in the critical 
temperature. 

AJVeb: Including higher bands brings additional states and 
hence increases the number of atoms in the saturated 
thermal cloud. This has the effect of decreasing the crit- 
ical temperature. 

Interestingly the dominant corrections at low temperatures 
(ATVle and AA^^) both lead to an increase in Tc, whereas the 
dominant correction at higher temperatures ( AA^eb) shifts Tc 
downwards. 



D. Numerical calculations of Tc 

While the analytic calculation provides a useful critical 
temperature estimate, the complexity of the spectrum in the 
combined harmonic -lattice potential necessitates a numerical 
solution. Here we discuss our procedure for calculating the 
critical temperature using the spectrum determined by full nu- 
merical diagonalization of ([T]i and give a simple numerical 
scheme that makes use of the piecewise density of states we 
have developed in Sees. IIIC l| and |IIC2| 



1. Full numerical calculation 




X 10 



Figure 4: Comparison of analytic and numerical critical temper- 
atures (a) Full numerical results for Tc (squares), analytic results 
Tco (dotted), Tci (dashed) and simple numerical result using piece- 
wise analytic density of states Tcjv (solid), (b) Energy scales com- 
pared to TcC). (c) ATV-corrections. Calculation parameters: isotropic 
harmonic trap with uo = Q.Q2^ujr, and lattice depth parameters 



Vj=&ERtor:3 = 1,2,3. 



From the results of our full diagonalization we deter- 
mine the one-dimensional energy spectrum {e^'} over 
a large energy range, typically including all states up 
to energy 25E]i above the ID ground state energy (as 
discussed in Sec. Ill Dl ). The thermal properties of the 
system are then calculated over a temperature range by 
iterating the chemical potential ji to find the desired total 
number of atoms, i.e. root-finding the expression/ (/i) — 

E., Jexp([e« + + e['^ - nykT) - 1}-^ n] 
for each T. From this calculation we hence evaluate the 
condensate population as a function of temperature, i.e, 
iVo(r) = {exp(E^ti - ^A/kT) - and de- 

termine the condensation temperature as that at which 
\{dNQ/dT)/NiA (i.e. the relative change in the ground state 
occupation) is maximised. 

2. Simple numerical calculation 

The critical temperature can also be estimated by perform- 
ing a simple numerical integral using the piecewise analytic 
density of states (flSl l under the saturated thermal cloud condi- 
tion {i.e. /i ^ Co) 



^(^)-£""^^ %(e-ef)/ll_l - (25) 

This result can then be numerically inverted to give a critical 
temperature estimate T^n — TcN {N) . The energy i^max ap- 
pearing in the integral has to be chosen such that Smax ^ kT, 
in which case the result will be independent of ii^max- 

This approach is significantly simpler than the full numer- 
ical calculation because it does not require a full numerical 
diagonalization. Indeed the information needed for g{e) can 
be obtained from results of the uniform lattice or tight-binding 
approximations, as discussed in Sec. IIIBI 

E. Comparison of analytic and numerical critical 
temperatures 

In Fig. Ua) we show analytic and numerical results for the 
critical temperature. To relate these parameters to those in ex- 
periment we note that for ^^Rb in a lattice with a — 425nm, 
the trap frequency corresponds to w 27r x 31 while the 
temperature scale is Eji/k w 152nK. These results show the 
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general behaviour we have observed over a wide parameter 
regime. Tco provides a useful critical temperature estimate, 
though is noticeably shifted relative to the full numerical re- 
sult. Including first order corrections Td provides a quanti- 
tatively much more accurately result, although its agreement 
with the full numerical result worsens for large N. Interest- 
ingly the simple numerical result Tcn outlined in Sec. IIII D 21 
provides an accurate description over the full range consid- 
ered. 

In Fig. |4fb) and (c) we explore the validity conditions for 
our derivation of the critical temperature. Note for the po- 
tential parameters used for the results in Fig. |4]we have that 
Ele - eo = Ele - wq = 0.123Er, ei, - wq = 

3.83Eii, and 62 — 60 — 6.7Er. The relative size of the param- 
eters {ELE-eo)/kTco, {ELE-wo)/kTco and kTco/{eij-wo) 
are shown in Fig. SJb). We require all of these parameters to 
be small for our analytic calculation to be valid. These re- 
sults show that for small N the critical temperature is suf- 
ficiently low that a first order treatment of the low energy 
spectrum is not appropriate (i.e. both (£^le — eo)/kTco and 
{Ele - wo)/kTco are large). 

At larger atom numbers (N) the term kTco / {eij ~wa) tends 
to grow reflecting the increased importance of excited band 
states. In Fig. HJc) we show the related values of AA'le, AA^^ 
and ANeb- At small N ( and hence small Tcq) the expansions 
we have used to obtainAA^LE and ATV^ are not valid. As N 
increases these contributions become less significant relative 
to N, however as AN^ scales like \i^(E[E^^^wo)Jkl^ it de- 
creases rather slowly with increasing T^q. Finally, the excited 
band contribution becomes gradually more significant with in- 
creasing number 



IV. GENERAL BEHAVIOUR OF CONDENSATION IN THE 
COMBINED POTENTIAL 

In figure |5] we show the results of our full numerical cal- 
culation (as discussed in Sec. IIII D 1 b for the critical temper- 
ature and condensate fraction over a wide parameter regime. 
For the case of uj — 0.025a'/; in Fig. 13 a) we see that as the 
lattice depth increases the critical temperature of the system 
decreases. While for the case of tD = O.ObcuR shown in Fig. 
ISjb) the critical temperature instead tends to increase with in- 
creasing lattice depth (for N sufficiently large). 

To understand these results we recall the critical tempera- 
ture for a harmonically trapped gas 

_fiLuf 2N \ ^ 

^^-^^--k[r{3)c{3)) • ^^^^ 

In comparison to our analytic result given in Eq. ( fTsl l. we note 
that the critical temperature for the combined potential scales 
with mean trap frequency and total atom number at higher 
powers, i.e. uj'^ and A^^/-^ respectively. Locating the trap fre- 
quency at which the critical temperatures for the harmonic and 
combined potentials are the same determines a critical mean 




Figure 5; Bose-Einstein condensation in a combined harmonic lat- 
tice potential, (a) Critical temperature as a function of lattice depth 
(all Vj = Vo) and total atom number for an isotropic harmonic trap 
with ijj = 0.025ci;_R. Isothermal levels spaced by 0.2ER/k shown 
as contour lines, (b) As for (a) but with lo = 0.05ujr. (c) Conden- 
sate fraction versus temperature for the same parameters as (a) and 
iV = 1 X 10^ atoms. 



trap frequency (uJc{N)y. 

cm _^ (airy 1 

COR 7r[ C(3) ) ^^'^ 
For ui > ujc the critical temperature is higher in the combined 
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potential than for the pure harmonic trap, whereas for (D < Wc 
the pure harmonic potential has a higher critical temperature. 
For N — 10^ atoms we find that uic ^ 0.049w/j, which is 
consistent with Figs. |5ja) and|5lb) which lie either side of this 
value. Since c^c is based on the simple critical temperature 
estimate (fTsT l. it will only be valid for cases where the critical 
temperature is not too high or low (as given by the validity 
conditions in Sec. IlllCb . 

In Fig. |3c) we show the condensate fraction versus temper- 
ature for a system of 10^ atoms in a combined potential with 
cZj = 0.025a'i^. As the lattice depth increases the critical tem- 
perature shifts downwards (as can also be discerned from Fig. 
lUa)), and the characteristic shape of the condensate fraction 
dependence on temperature, A^o ~ [1 — (T/Tc)"], changes 
from Q!~3tOQ;~3/2. These predicted features should be 
verifiable by current experiments. 

V. CONCLUSION 

We have performed a comprehensive study of the critical 
temperature for an ideal Bose gas in a combined harmonic 
lattice potential. We have described distinctive regions of 
the spectrum and have shown that a simple piecewise density 



of states provides an accurate characterization of this system 
for lattice depths greater than about iEj^. We have devel- 
oped an analytic expression for the critical temperature in the 
combined potential. The corrections to this result are typi- 
cally significant, and we have shown that including them pro- 
vides a useful estimate for the critical temperature obtained 
by a full numerical calculation. Additionally, we give a sim- 
ple numerical procedure based on piecewise density of states 
that provides an accurate prediction for the critical tempera- 
ture. Finally we have presented results over a wide parame- 
ter regime appropriate to current experiments and have shown 
that the critical temperature in the combined potential can be 
increased or decreased relative to that of the pure harmonic 
trap. 
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